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Abstract: The focus of this article is on entropy and Markov processes. We study the 
properties of functional which are invariant with respect to monotonic transformations and 
analyze two invariant "additivity" properties: (i) existence of a monotonic transformation 
which makes the functional additive with respect to the joining of independent systems 
and (ii) existence of a monotonic transformation which makes the functional additive with 
respect to the partitioning of the space of states. All Lyapunov functionals for Markov 
chains which have properties (i) and (ii) are derived. We describe the most general ordering 
of the distribution space, with respect to which all continuous-time Markov processes are 
monotonic (the Markov order). The solution differs significantly from the ordering given by 
the inequality of entropy growth. For inference, this approach results in a convex compact 
set of conditionally "most random" distributions. 

Keywords: Markov process; Lyapunov function; entropy functionals; attainable region; 
MaxEnt; inference 



1. Introduction 

1.1. A Bit of History: Classical Entropy 

Two functions, energy and entropy, rule the Universe. 
In 1865 R. Clausius formulated two main laws [1]: 
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1 . The energy of the Universe is constant. 

2. The entropy of the Universe tends to a maximum. 

The universe is isolated. For non-isolated systems energy and entropy can enter and leave, the change 
in energy is equal to its income minus its outcome, and the change in entropy is equal to entropy 
production inside the system plus its income minus outcome. The entropy production is always positive. 

Entropy was born as a daughter of energy. If a body gets heat AQ at the temperature T then for this 
body dS = AQ/T. The total entropy is the sum of entropies of all bodies. Heat goes from hot to cold 
bodies, and the total change of entropy is always positive. 

Ten years later J.W. Gibbs [2] developed a general theory of equilibrium of complex media based on 
the entropy maximum: the equilibrium is the point of the conditional entropy maximum under given 
values of conserved quantities. The entropy maximum principle was applied to many physical and 
chemical problems. At the same time J.W. Gibbs mentioned that entropy maximizers under a given 
energy are energy minimizers under a given entropy. 

The classical expression J plnp became famous in 1872 when L. Boltzmann proved his 
if-theorem [3]: the function 

H = J f{x,v)\nf{x,v)dxdv 

decreases in time for isolated gas which satisfies the Boltzmann equation (here f{x,v) is the distribution 
density of particles in phase space, x is the position of a particle, v is velocity). The statistical entropy 
was born: S = —kH. This was the one-particle entropy of a many -particle system (gas). 

In 1902, J.W. Gibbs published a book "Elementary principles in statistical dynamics" [4]. He 
considered ensembles in the many-particle phase space with probability density p{pi,qi, ■ ■ ■Pn,Qn), 
where pi, Qi are the momentum and coordinate of the ith particle. For this distribution, 

S = -k J p(pi,gi,...p„,g„)ln(p(pi,gi,...p„,g„))dgi...dg„dpi...dp„ (1) 

Gibbs introduced the canonical distribution that provides the entropy maximum for a given expectation 
of energy and gave rise to the entropy maximum principle (MaxEnt). 

The Boltzmann period of history was carefully studied [5]. The difference between the Boltzmann 
entropy which is defined for coarse-grained distribution and increases in time due to gas dynamics, and 
the Gibbs entropy, which is constant due to dynamics, was analyzed by many authors [6,7]. Recently, 
the idea of two functions, energy and entropy which rule the Universe was implemented as a basis of 
two-generator formalism of nonequilibrium thermodynamics [8,9]. 

In information theory, R.V.L. Hartley (1928) [10] introduced a logarithmic measure of information 
in electronic communication in order "to eliminate the psychological factors involved and to establish 
a measure of information in terms of purely physical quantities". He defined information in a text of 
length n in alphabet of s symbols as H = n log s. 

In 1948, C.E. Shannon [11] generalized the Hartley approach and developed "a mathematical theory 
of communication", that is information theory. He measured information, choice and uncertainty by 
the entropy: 

n 

S = -"^Pilogpi (2) 

i=l 
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Here, Pi are the probabilities of a full set of n events (Yl^^iPi = 1). The quantity 5* is used to 
measure of how much "choice" is involved in the selection of the event or of how uncertain we are 
of the outcome. Shannon mentioned that this quantity form will be recognized as that of entropy, as 
defined in certain formulations of statistical mechanics. The classical entropy (1), (2) was called the 
Boltzmann-Gibbs-Shannon entropy (BGS entropy). (In 1948, Shannon used the concave function (2), 
but under the same notation H as for the Boltzmann convex function. Here we use H for the convex 
_fr-function, and S for the concave entropy.) 

In 1951, S. KuUback and R.A. Leibler [12] supplemented the BGS entropy by the relative BGS 
entropy, or the KuUback-Leibler divergence between the current distribution P and some "base" (or 
"reference") distribution Q: 

D^^{P\\Q) = Y,P^^og^- (3) 

i 

The KuUback-Leibler divergence is always non-negative _Dkl(-P||Q) > (the Gibbs inequality). It 
is not widely known that this "distance" has a very clear physical interpretation. This function has been 
well known in physical thermodynamics since 19th century under different name. If Q is an equilibrium 
distribution at the same temperature as P has, then 

where F is free energy and T is thermodynamic temperature. In physics, F = U — TS, where 
physical entropy S includes an additional multiplier k, the Boltzmann constant. The thermodynamic 
potential —F/T has its own name, Massieu function. Let us demonstrate this interpretation of the 
KuUback-Leibler divergence. The equilibrium distribution Q provides the conditional entropy (2) 
maximum under a given expectation of energy ^^Uiqi = U and the normalization condition Xli 5'« = 1- 
With the Lagrange multipliers fiu and /io we get the equilibrium Boltzmann distribution: 

expj-fiuUj) 

Qi = exp(-/io - iJ^uUi) = ^ 7 r (5) 

The Lagrange multiplier fiu is in physics (by definition) 1 / kT, so S{Q) = + hence, j^o = 
For the KuUback-Leibler divergence this formula gives (4). 

After the classical work of Zeldovich (1938, reprinted in 1996 [13]), the expression for free energy in 
the "KuUback-Leibler form" 

Ci 



where Cj is concentration and c*(T) is the equilibrium concentration of the ith component, is recognized 
as a useful instrument for the analysis of kinetic equations (especially in chemical kinetics [14,15]). 

Each given positive distribution Q could be represented as an equilibrium Boltzmann distribution for 
given T > if we take Ui = —kT log + uq for an arbitrary constant level uq. If we change the order 
of arguments in the KuUback-Leibler divergence then we get the relative Burg entropy [16,17]. It has 
a much more exotic physical interpretation: for a current distribution P we can define the "auxiliary 
energy" functional Up for which P is the equilibrium distribution under a given temperature T. We can 
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calculate the auxiliary free energy of any distribution Q and this auxiliary energy functional: Fp{Q). (Up 
to an additive constant, for P = P* this Fp(Q) turns into the classical free energy, Fp{Q) = F(Q).) In 
particular, we can calculate the auxiliary free energy of the physical equilibrium, Fp(P*). The relative 
Burg entropy is 

This functional should also decrease in any Markov process with given equilibrium P*. 

Information theory developed by Shannon and his successors focused on entropy as a measure of 
uncertainty of subjective choice. This understanding of entropy was returned from information theory to 
statistical mechanics by E.T. Jaynes as a basis of "subjective" statistical mechanics [18,19]. He followed 
Wigner's idea "entropy is an antropocentric concept". The entropy maximum approach was declared 
as a minimization of the subjective uncertainty. This approach gave rise to a MaxEnt "anarchism". It 
is based on a methodological hypothesis that everything unknown could be estimated by the principle 
of the entropy maximum under the condition of fixed known quantities. At this point the classicism 
in entropy development changed to a sort of scientific modernism. The art of model fitting based on 
entropy maximization was developed [20] . The principle of the entropy maximum was applied to plenty 
of problems: from many physical problems [21], chemical kinetics and process engineering [15] to 
econometrics [22,23] and psychology [24]. Many new entropies were invented and now one has rich 
choice of entropies for fitting needs [25]. The most celebrated of them are the Renyi entropy [26], the 
Burg entropy [16,17], the Tsallis entropy [27,28] and the Cressie-Read family [29,30]. The nonlinear 
generalized averaging operations and generalized entropy maximization procedures were suggested [31]. 

Following this impressive stream of works we understand the MaxEnt approach as conditional 
maximization of entropy for the evaluation of the probability distribution when our information is partial 
and incomplete. The entropy function may be the classical BGS entropy or any function from the rich 
family of non-classical entropies. This rich choice causes a new problem: which entropy is better for a 
given class of applications? 

The MaxEnt "anarchism" was criticized many times as a "senseless fitting". Arguments pro and 
contra the MaxEnt approach with non-classical entropies (mostly the Tsallis entropy [27]) were collected 
by Cho [32]. This sometimes "messy and confusing situation regarding entropy-related studies has 
provided opportunities for us: clearly there are still many very interesting studies to pursue" [33]. 

1.2. Key Points 

In this paper we do not pretend to invent new entropies. (There appear new functions as limiting cases 
of the known entropy families, but this is not our main goal). Entropy is understood in this paper as a 
measure of uncertainty which increases in Markov processes. In our paper we consider a Markov process 
as a semigroup on the space of positive probability distributions. The state space is finite. Generalizations 
to compact state spaces are simple. We analyze existent relative entropies (divergences) using several 
simple ideas: 



1. In Markov processes probability distributions P{t) monotonically approach equilibrium P*: 
divergence D{P{t)\\P*) monotonically decrease in time. 
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2. In most applications, conditional minimizers and maximizers of entropies and divergences are 
used, but the values are not. This means that the system of level sets is more important than the 
functions' values. Hence, most of the important properties are invariant with respect to monotonic 
transformations of entropy scale. 

3. The system of level sets should be the same as for additive functions: after some rescaling the 
divergences of interest should be additive with respect to the joining of statistically independent 



4. The system of level sets should after some rescaling the divergences of interest should have the 
form of a sum (or integral) over states J2i fiPi^Pi)^ where the function / is the same for all 
states. In information theory, divergences of such form are called separable, in physics the term 
trace-form functions is used 

The first requirement means that if a distribution becomes more random then it becomes closer 
to equilibrium (Markov process decreases the information excess over equilibrium). For example, 
classical entropy increases in all Markov processes with uniform equilibrium distributions. This is 
why we may say that the distribution with higher entropy is more random, and why we use entropy 
conditional maximization for the evaluation of the probability distribution when our information is partial 
and incomplete. 

It is worth to mention that some of the popular Bregman divergences, for example, the squared 
Euclidean distance or the Itakura-Saito divergence, do not satisfy the first requirement (see Section 4.3). 

The second idea is just a very general methodological thesis: to evaluate an instrument (a divergence) 
we have to look how it works (produces conditional minimizers and maximizers). The properties 
of the instrument which are not related to its work are not important. The requirements number 
three and four are, in their essence, conditions for separation of variables in conditional minimization 
(maximization) problems. 

The number three allows to separate variables if the system consists of independent subsystems, the 
number four relates to separation of variables for partitions of the space of probability distributions. 

Amongst a rich world of relative entropies and divergences, only two families meet these 
requirements. Both were proposed in 1984. The Cressie-Read (CR) family [29,30]: 



and the convex combination of the Burg and Shannon relative entropies proposed in [34] and further 
analyzed in [35,36]: 



When A —7- the CR divergence tends to the KL divergence (the relative Shannon entropy) and when 
A — 7- —1 it turns into the Burg relative entropy. The Tsallis entropy was introduced four years later 
[27] and became very popular in thermostatistics (there are thousands of works that use or analyze this 



systems. 




A g] — oo, oo[ 
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entropy [37]). The Tsallis entropy coincides (up to a constant multiplier A + 1) with the CR entropy for 
A g] — 1, oo[ and there is no need to study it separately (see discussion in Section 2.2). 

A new problem arose: which entropy is better for a specific problem? Many authors compare 
performance of different entropies and metrics for various problems (see, for example, [39,40]). In 
any case study it may be possible to choose "the best" entropy but in general we have no sufficient 
reasons for such a choice. We propose a possible way to avoid the choice of the best entropy. 

Let us return to the idea: the distribution Q is more random than P if there exists a continuous-time 
Markov process (with given equilibrium distribution P*) that transforms P into Q. We say in this case 
that P and Q are connected by the Markov preorder with equilibrium P* and use notation P y-^p, Q. 
The Markov order yp* is the transitive closure of the Markov preorder. 

If a priori information gives us a set of possible distributions W then the conditionally "maximally 
random distributions" (the "distributions without additional information", the "most indefinite 
distributions" in W) should be the minimal elements in W with respect to Markov order. If a Markov 
process (with equilibrium P*) starts at such a minimal element P then it cannot produce any other 
distribution from W because all distributions which are more random that P are situated outside W. In 
this approach, the maximally random distributions under given a priori information may be not unique. 
Such distributions form a set which plays the same role as the standard MaxEnt distribution. For the 
moment based a priori information the set W is an intersection of a linear manifold with the simplex 
of probability distributions, the set of minimal elements in W is also polyhedron and its description is 
available in explicit form. In low-dimensional case it is much simpler to construct this polyhedron than 
to find the MaxEnt distributions for most of specific entropies. 

1.3. Structure of the Paper 

The paper is organized as follows. In Section 2 we describe the known non-classical divergences 
(relative entropies) which are the Lyapunov functions for the Markov processes. We discuss the general 
construction and the most popular families of such functions. We pay special attention to the situations, 
when different divergences define the same order on distributions and provide the same solutions of the 
MaxEnt problems (Section 2.2). In two short technical Sections 2.3 and 2.4 we present normalization 
and symmetrization of divergences (similar discussion of these operations was published very recently 
[38]. 

The divergence between the current distribution and equilibrium should decrease due to Markov 
processes. Moreover, divergence between any two distributions should also decrease (the generalized 
data processing Lemma, Section 3). 

Definition of entropy by its properties is discussed in Section 4. Various approaches to this definition 
were developed for the BGS entropy by Shannon [11], [41] and by other authors for the Renyi entropy 
[43,44], the Tsallis entropy [42], the CR entropy and the convex combination of the BGS and Burg 
entropies [46]. Csiszar [45] axiomatically characterized the class of Csiszar-Morimoto divergences 
(formula (6) below). Another characterization of this class was proved in [46] (see Lemma 1, Section 4.3 
below). 

From the celebrated properties of entropy [47] we selected the following three: 



Entropy 2010, 12 



1151 



1. Entropy should be a Lyapunov function for continuous -time Markov processes; 

2. Entropy is additive with respect to the joining of independent systems; 

3. Entropy is additive with respect to the partitioning of the space of states {i.e., has the trace-form). 

To solve the MaxEnt problem we have to find the maximizers of entropy (minimizers of the relative 
entropy) under given conditions. For this purpose, we have to know the sublevel sets of entropy, 
but not its values. We consider entropies with the same system of sublevel sets as equivalent ones 
(Section 2.2). From this point of view, all important properties have to be invariant with respect to 
monotonic transformations of the entropy scale. Two last properties from the list have to be substituted 
by the following: 

2'. There exists a monotonic transformation which makes entropy additive with respect to the joining 
of independent systems (Section 4.1); 

3'. There exists a monotonic transformation which makes entropy additive with respect to the 
partitioning of the space of states (Section 4.2). 

These properties imply specific separation of variables for the entropy maximization problems 
(Sections 4.1 and 4.2). Several "No More Entropies" Theorems are proven in Section 4.3: if an entropy 
has properties 1, 2' and 3' then it belongs to one of the following one-parametric families: to the 
Cressie-Read family, or to a convex combination of the classical BGS entropy and the Burg entropy 
(may be, after a monotonic transformation of scale). 

It seems very natural to consider divergences as orders on distribution spaces (Section 5.1), the 
sublevel sets are the lower cones of this orders. For several functions, Hi{P), . . . ,Hk{P) the sets 
{Q I Hi{P) > Hi{Q) for all i} give the simple generalization of the sublevel sets. In Section 5 we 
discuss the more general orders in which continuous time Markov processes are monotone, define the 
Markov order and fully characterize the local Markov order. The Markov chains with detailed balance 
define the Markov order for general Markov chains (Section 5.2). It is surprising that there is no necessity 
to consider other Markov chains for the order characterization (Section 5.2) because no reversibility is 
assumed in this analysis. 

In Section 6. 1 we demonstrate how is it possible to use the Markov order to reduce the uncertainty in 
the standard settings when a priori information is given about values of some moments. Approaches to 
construction of the most random distributions are presented in Section 6.2. 

Various approaches for the definition of the reference distribution (or the generalized canonical 
distribution) are compared in Section 7. 

In Conclusion we briefly discuss the main results. 
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2. Non- Classical Entropies 

2.1. The Most Popular Divergences 

Csiszar-Morimoto Functions Hh 

During the time of modernism plenty of new entropies were proposed. Esteban and Morales [25] 
attempted to systemize many of them in an impressive table. Nevertheless, there are relatively few 
entropies in use now. Most of the relative entropies have the form proposed independently in 1963 by 
I. Csiszar [49] and T. Morimoto [48]: 



where h{x) is a convex function defined on the open (x > 0) or closed x > semi-axis. We use here 
notation Hh{P\\P*) to stress the dependence of both on pi and p*. 

These relative entropies are the Lyapunov functions for all Markov chains with equilibrium 
P* = (p*). Moreover, they have the relative entropy contraction property given by the generalized 
data processing lemma (Section 3.2 below). 

For h(x) = xlogx this function coincides with the KuUback-Leibler divergence from the current 
distribution pi to the equilibrium p*. Some practically important functions h have singularities at 
0. For example, if we take h{x) = — logx, then the correspondent Hh is the relative Burg entropy 



Required Properties of the Function h{x) 

Formally, h{x) is an extended real-valued proper convex function on the closed positive real half-line 
[0, oo[. An extended real-valued function can take real values and infinite values ±oo. A proper function 
has at least one finite value. An extended real valued function on a convex set U is called convex if 
its epigraph 



is a convex set [50]. For a proper function this definition is equivalent to the Jensen inequality 



It is assumed that the function h{x) takes finite values on the open positive real half-line ]0, oo[ but 
the value at point x = may be infinite. For example, functions h{x) = — log x or h(x) = {a > 0) 
are allowed. A convex function h{x) with finite values on the open positive real half-line is continuous 
on ]0, oo[ but may have a discontinuity at x = 0. For example, the step function, h{x) = if x = and 
h{x) = — 1 if X > 0, may be used. 

A convex function is differentiable almost everywhere. Derivative of h{x) is a monotonic function 
which has left and right limits at each point x > 0. An inequality holds: h'{x){y — x) < h(y) — h{x) 
(Jensen's inequality in the differential form). It is valid also for left and right limits of h! at any point 
X > 0. 




(6) 



Hh = - J2iPi logfe/p*) ^ oo for Pi 0. 



epi{h) = {(x,y) I X > 0,y > h{x)} 



h{ax + (1 — a)y) < ah{x) + (1 — a)h{y) for all x,y e U, a e [0, 1] 
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Not everywhere differentiable functions h{x) are often used, for example, h{x) = |x — 1|. 
Nevertheless, it is convenient to consider the twice differentiable on ]0, oo[ functions h(x) and to produce 
a non-smooth h(x) (if necessary) as a limit of smooth convex functions. We use widely this possibility. 

The Most Popular Divergences Hh{P\\P*) 

1. Let h{x) be the step function, h{x) = if x = and h{x) = — 1 if x > 0. In this case, 

Hj,iP\\P*) = - 5^ 1 (7) 

i, pi>0 

The quantity —Hh is the number of non-zero probabilities Pi and does not depend on P*. 
Sometimes it is called the Hartley entropy. 

2. h = \x — 1\, 

Hn{P\\P*) = J2\p.-p*\ 

i 

this is the /i -distance between P and P*. 

3. h = xlnx. 



Hh{P\\P*) = J^P.ln = D^^{P\\P*) 



(8) 



this is the usual KuUback-Leibler divergence or the relative BGS entropy; 
4. h = — hix. 



Hh{P\\P*) = (§) = Dkl{P* 



\\P) (9) 



this is the relative Burg entropy. It is obvious that this is again the KuUback-Leibler divergence, 
but for another order of arguments. 

5. Convex combinations of h = xlnx and h = — Inx also produces a remarkable family of 
divergences: h = /3x\nx — {1 — /3)\nx (/3 E [0, 1]), 

Hh{P\\P*) = PDKLiP\\P*) + (1 - f3)DKL{P*\\P) (10) 

The convex combination of divergences becomes a symmetric functional of (P, P*) for /3 = 1/2. 
There exists a special name for this case, "Jeffreys' entropy". 



6. h 



Hh{P\\P*) = IY1 



(11) 



Pi 

This is the quadratic term in the Taylor expansion of the relative Botzmann-Gibbs-Shannon 
entropy, DKLiP\\P*), near equilibrium. Sometimes, this quadratic form is called the Fisher 
entropy. 
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1 h - ^(^^"1) 

A(A+1) ' 



Hh{P\\P*) 



A(A 



2^ P. 



^1 -1 

Pi 



(12) 



This is the CR family of power divergences [29,30]. For this family we use notation i^cR a- If 
A then Hck\ -Dkl(-P||-P*), this is the classical BGS relative entropy; if A -1 then 
HcK A -Dkl(-P*||-P), this is the relative Burg entropy. 

8. For the CR family in the limits A — j- ±oo only the maximal terms "survive". Exactly as we get the 
limit /°° of F norms for p — )■ oo, we can use the root (A(A + l)-ffcR a)^^''^' for A — )■ ±oo and write 
in these limits the divergences: 



i^CRoo(P|m=max|p|-l 



^^CR-oo(P||P*) = max<{^[.-l 



(13) 



(14) 



The existence of two limiting divergences ifcR ±oo seems very natural: there may be two types 
of extremely non-equilibrium states: with a high excess of current probability pi above p* and, 
inversely, with an extremely small current probability Pi with respect to p*. 

9. The Tsallis relative entropy [27] corresponds to the choice h = -^^^f^p, a > 0. 

H,{P\\P*) = ^J2p^ 



a 

For this family we use notation i^xs a- 
Renyi Entropy 

The Renyi entropy of order a > is defined [26] as 

It is a concave function, and 



(15) 



(p)^^iog(|:.r) 



(16) 



Hno.{P)^S{P) 

when a —J- 1, where S{P) is the Shannon entropy. 

When a — )■ oo, the Renyi entropy has a limit Hoo{X) = — logmaxj=i^...„pj, which has a special 
name "Min-entropy". 

It is easy to get the expression for a relative Renyi entropy H^i a(P||P*) from the requirement that it 
should be a Lyapunov function for any Markov chain with equilibrium P*: 



a-r 
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For the Min-entropy, the correspondent divergence (the relative Min-entropy) is 

i7,o(P|m = log max (^] 

It is obvious from (22) below that maxj=i is a Lyapunov function for any Markov chain with 

equilibrium P* , hence, the relative Min-entropy is also the Lyapunov functional. 

2.2. Entropy Level Sets 

A level set of a real-valued function / of is a set of the form : 

{x I /(x) = c} 

where c is a constant (the "level"). It is the set where the function takes on a given constant value. A 
sublevel set of / is a set of the form 

{x I f{x) < c} 

A superlevel set of / is given by the inequality with reverse sign: 

{x I /(x) > c} 

The intersection of the sublevel and the superlevel sets for a given value c is the level set for this level. 

In many applications, we do not need the entropy values, but rather the order of these values on the 
line. For any two distributions P, Q we have to compare which one is closer to equilibrium P*, i.e., to 
answer the question: which of the following relations is true: H{P\\P*) > H{Q\\P*), H{P\\P*) = 
H{Q\\P*) or H{P\\P*) < H{Q\\P*)1 To solve the MaxEnt problem we have to find the maximizers of 
entropy (or, in more general settings, the minimizers of the relative entropy) under given conditions. For 
this purpose, we have to know the sublevel sets, but not the values. All the MaxEnt approach does not 
need the values of the entropy but the sublevel sets are necessary. 

Let us consider two functions, and ^ on a set U. For any V C U we can study conditional 
minimization problems 0(x) min and 'ip{x) — )■ min, x E V. The sets of minimizers for these two 
problems coincide for any V C U if and only if the functions (p and have the same sets of sublevel 
sets. It should be stressed that here just the sets of sublevel sets have to coincide without any relation to 
values of level. 

Let us compare the level sets for the Renyi, the Cressie-Read and the Tsallis families of divergences 
(for a — 1 = X and for all values of a). The values of these functions are different, but the level sets are 
the same (outside the Burg limit, where a 0): for a 7^ 0, 1 

H^^{P\\P*) = ^—\nc- HcKa-iiP\\P*) = ^^-^{c-iy, HrsaiP\\P*) = ^-^ic-l) (17) 
a — 1 a{a — I) a — 1 

where c = Y.iPiiPi/PiT~^- 
Beyond points a = 0, 1 

HcKa-iiP\\P*) = —^-—exp{{a-l)Hna{P\\P*)) = -Ht,^{P\\P*) 
a[a — 1 a 
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For a — > 1 all these divergences turn into the Shannon relative entropy. Hence, if a 7^ then for any 
P, P*, Q, Q* the following equalities A, B, C are equivalent, A^B'^C: 



This equivalence means that we can select any of these three divergences as a basic function and 
consider the others as functions of this basic one. 

For any a > and \ = a + 1 the Renyi, the Cressie-Read and the Tsallis divergences have the 
same family of sublevel sets. Hence, they give the same maximizers to the conditional relative entropy 
minimization problems and there is no difference which entropy to use. 

The CR family has a more convenient normalization factor 1/A(A + 1) which gives a proper convexity 
for all powers, both positive and negative, and provides a sensible Burg limit for A — j- —1 (in contrary, 
when a —J- both the Renyi and Tsallis entropies tend to 0). 

When a < then for the Tsallis entropy function h = loses convexity, whereas for the 

Cressie-Read family convexity persists for all values of A. The Renyi entropy also loses convexity for 
a < 0. Neither the Tsallis, nor the Renyi entropy were invented for use with negative a. 

There may be a reason: for a < the function x" is defined for a: > only and has a singularity at x = 
0. If we assume that the divergence should exist for all non-negative distributions, then the cases a < 
should be excluded. Nevertheless, the Burg entropy which is singular at zeros is practically important and 
has various attractive properties. The Jeffreys' entropy (the symmetrized KuUback-Leibler divergence) 
is also singular at zero, but has many important properties. We can conclude at this point that it is not 
obvious that we have to exclude any singularity at zero probability. It may be useful to consider positive 
probabilities instead ("nature abhors a vacuum"). 

Finally, for the MaxEnt approach (conditional minimization of the relative entropy), the Renyi and 
the Tsallis families of divergences (a > 0) are particular cases of the Cressie-Read family because they 
give the same minimizers. For a < the Renyi and the Tsallis relative entropies lose their convexity, 
while the Cressie-Read family remains convex for A < —1 too. 

2.3. Minima and normalization 

For a given P*, the function Hh achieves its minimum on the hyperplane XliP* = ZliK =const at 
equilibrium p*, because at this point 



The transformation h{x) — h{x) + ax + 6 just shifts Hh by constant value: Hh — > Hh + a YliiPi + b = 
Hh + a + h. Therefore, we can always assume that the function h{x) achieves its minimal value at point 
X = 1, and this value is zero. For this purpose, one should just transform h: 



A. H^^{P\\P*) = H^^{Q\\Q*) 

B. HcKa+l{P\\P*) = HcKa+imQl 

C. i/Tsa(P|m=^Tsa(Q||g*) 



(18) 




h{x) ■.= h{x)-h{l)-h'{l){x-l) 



(19) 
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This normalization transforms xlnx into xlnx — (x — 1), — Inx into — Inx + (x — 1), and into 
x"" — 1 — a{x — \). After normalization Hh{P\\P*) > 0. If the normalized h{x) is strictly positive 
outside point x = 1 {h{x) > if x 7^ 1) then Hh{P\\P*) = if and only if P = P* {i.e., in equilibrium). 

The normalized version of any divergence Hh{P\\P*) could be produced by the normalization 
transformation h{x) := h{x) — h{l) — h'(l){x — 1) and does not need separate discussion. 

2.4. Symmetrization 

Another technical issue is symmetry of a divergence. If h{x) = xlnx then both Hh{P\\P*) (the 
KL divergence) and Hh{P*\\P) (the relative Burg entropy) are the Lyapunov functions for the Markov 
chains, and Hh{P*\\P) = Hg{P\\P*) with g{x) = — Inx. Analogously, for any h{x) we can write 
iJ^(P*||P) =iJ,(P||P*) with 

g{x) = xh (-] (20) 



X ^ 

If h{x) is convex on R+ then g{x) is convex on R+ too because 

g"{x) = -h"^^ 



x^ \x 



The transformation (20) is an involution: 



^9 [ - ] = K^) 



X 



The fixed points of this involution are such functions h{x) that Hh{P\\P*) is symmetric with respect 
to transpositions of P and P*. There are many such h{x). An example of symmetric Hh{P\\P*) gives 
the choice h{x) = —\fx\ 

u,iP\\p^) = -Y,^, 

i 

After normalization, we get 

h{x):=\{^x-\f- H,{P\\P*) = lj2(y^^- 



Essentially (up to a constant addition and multiplier) this function coincides with a member of the CR 
family, Hq^i _i (12), and with one of the Tsallis relative entropies Hr^,^ 1 (15). The involution (20) is a 
linear operator, hence, for any convex h{x) we can produce its symmetrization: 



]^{h{x)+g{x)) = ]^U{x)+xhr- 



For example, if h{x) = xlogx then hsymix) = ^{xlogx — logx); if h{x) = x" then hsymix) 



i(x" + xi-") 
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3. Entropy Production and Relative Entropy Contraction 

3.1. Lyapunov Functionals for Markov Chains 

Let us consider continuous time Markov chains with positive equilibrium probabilities p*. The 
dynamics of the probability distribution pi satisfy the Master equation (the Kolmogorov equation): 

= ^{(lijPi - (21) 

where coefficients qij (i ^ j) are non-negative. For chains with a positive equilibrium distribution p* 
another equivalent form is convenient: 

dpi ^ * / Pj Pi \ 
where p* and q^j are connected by identity 

The time derivative of the Csiszar-Morimoto function Hh{p) (6) due to the Master equation is 



dHh{P\\P* 



dt 



5Z 1'iP*3 



-M^ +/^'^ p-^ 
p*ij \p]j \PiJ \Pj p*i 



< (24) 



To prove this formula, it is worth to mention that for any n numbers hi, j j^iQijP'jihi — hj) = 
0. The last inequality holds because of the convexity of h(x): h'(x){y — x) < h{y) — h{x) 
(Jensen's inequality). 
Liversely, if 

h{x)-h{y) + h'{y){x-y)<0 (25) 

for all positive x, y then h{x) is convex on R+. Therefore, if for some function h{x) Hh{p) is the 
Lyapunov function for all the Markov chains with equilibrium P* then h{x) is convex on R+. 

The Lyapunov functionals Hh do not depend on the kinetic coefficients q^j directly. They depend 
on the equilibrium distribution p* which satisfies the identity (23). This independence of the kinetic 
coefficients is the universality property. 

3.2. "Lyapunov Divergences" for Discrete Time Markov Chains 

The Csiszar-Morimoto functions (6) are also Lyapunov functions for discrete time Markov chains. 
Moreover, they can serve as a "Lyapunov distances" [51] between distributions which decreases due to 
time evolution (and not only the divergence between the current distribution and equilibrium). In more 
detail, let A = (aij) be a stochastic matrix in columns: 

O'ij > 0, ttij = 1 for all j 
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The ergodicity contraction coefficient for A is a number a{A) [52,53]: 

a(^) = \ max | ^ \aij - a^j \ | 

< a{A) < 1. 

Let us consider in this subsection the normalized Csiszar-Morimoto divergences Hh(P\\Q) (19): 
h{l) = 0, h{x) > 0. 

Theorem about relative entropy contraction. (The generalized data processing Lemma.) For each 
two probability positive distributions P, Q the divergence Hh{P\\Q) decreases under action of stochastic 
matrix A [54,55]: 

Hh{AP\\AQ) < a{A)Hh{P\\Q) (26) 



The generalizations of this theorem for general Markov kernels seen as operators on spaces 
of probability measures was given by [56]. The shift in time for continuous-time Markov chain 
is a column-stochastic matrix, hence, this contraction theorem is also valid for continuous-time 
Markov chains. 

The question about a converse theorem arises immediately. Let the contraction inequality hold for 
two pairs of positive distributions (P, Q) and (f/, V) and for all Hh'. 

HhiUWV) < Hh{P\\Q) (27) 

Could we expect that there exists such a stochastic matrix A that U = AP and V = API The answer 
is positive: 

The converse generalized data processing lemma. Let the contraction inequality (27) hold for two 
pairs of positive distributions (P, Q) and {U, V) and for all normalized H^. Then there exists such a 
column-stochastic matrix A that U = AP and V = AQ [54]. 

This means that for the system of inequalities (27) (for all normalized convex functions /i on ]0, oo[) 
is necessary and sufficient for existence of a (discrete time) Markov process which transform the pair 
of positive distributions (P, Q) in (U, V). It is easy to show that for continuous-time Markov chains 
this theorem is not valid: the attainable regions for them are strictly smaller than the set given by 
inequalities (27) and could be even non-convex (see [60] and Section 8.1 below). 

4. Definition of Entropy by its Properties 

4. 1 . Additivity Property 

The additivity property with respect to joining of subsystems is crucial both for the classical 
thermodynamics and for the information theory. 

Let us consider a system which is result of joining of two subsystems. A state of the system is an 
ordered pair of the states of the subsystems and the space of states of the system is the Cartesian product 
of the subsystems spaces of state. For systems with finite number of states this means that if the states 
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of subsystems are enumerated by indexes j and k then the states of the system are enumerated by pairs 
jk. The probability distribution for the whole system is pj^, and for the subsystems the probability 
distributions are the marginal distributions qj = '^i^Pjk, rk = J2jPjk- 

The additive functions of state are defined for each state of the subsystems and for a state of the whole 
system they are sums of these subsystem values: 

Uj^. = Vj + Wk 

where Vj and Wk are functions of the subsystems state. 

In classical thermodynamics such functions are called the extensive quantities. For expected values 
of additive quantities the similar additivity condition holds: 

UjkPjk = + Wk)Pik = ^ vjqj + ^ WkTk (28) 

j,k j,k j k 

Let us consider these expected values as functional of the probability distributions: u{{pjk}), v{{qj}) 
and w{{rk}). Then the additivity property for the expected values reads: 

<{P,k}) = v{{qj}) + w{{r,}) (29) 

where qj and the are the marginal distributions. 

Such a linear additivity property is impossible for non-linear entropy functionals, but under some 
independence conditions the entropy can behave as an extensive variable. 

Let P be a product of marginal distributions. This means that the subsystems are statistically 
independent: pjk = qjVk- Assume also that the distribution P* is also a product of marginal 
distributions p*;, = g^r^. Then some entropies reveal the additivity property with respect to joining of 
independent systems. 

1. TheBGS relative entropy Dki.{P\\P*) = DkUQUQ*) + DkURWR*)- 

2. The Burg entropy DKLiP*\\P) = P'kl(QIIQ) + DKLiR*\\R) ■ It is obvious that a convex 
combination of the Shannon and Burg entropies has the same additivity property. 

3. The Renyi entropy a{P\\P*) = Hr a{Q\\Q*) + ^^r aiR\\R*)- For a ^ oo the Min-entropy 
also inherites this property. 

This property implies the separation of variables for the entropy maximization problems. Let 
functionals u^{{pjk}), ■ ■ .u"^{{pjk}) be additive (28) (29) and let the relative entropy H{P\\P*) be 
additive with respect to joining of independent systems. Then the solution to the problem 

H{P\\P*) min 

subject to conditions 

u\P) = Ui (i = l,...m) 
is = g™'°r™'°, where g™™, r^^^ are solutions of partial problems: 

H{Q\\Q*) min 
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subject to the conditions 

v\Q) = V, (z = l,...m) 

and 

H{R\\R*) min 

subject to the conditions 

w\Q) = Wi (2 = l,...m) 

for some redistribution of the additive functionals values Ui = Vi + Wi. 

Let us call this property the separation of variables for independent subsystems. 

Neither the CR, nor the Tsallis divergences families have the additivity property. It is proven [46] that 
a function Hh has the additivity property if and only if it is a convex combination of the Shannon and 
Burg entropies. See also Theorem 3 in Appendix. 

Nevertheless, they have the property of separation of variables for independent subsystems because 
of the coincidence of the level sets with the additive function, the Renyi entropy (for all a > 0). 

The Tsallis entropy family has absolutely the same property of separation of variables as the Renyi 
entropy. To extend this property of the Renyi Tsallis entropies for negative a, we have to change there 
min to max. 

For the CR family the result sounds even better: because of better normalization, the separation of 
variables is valid for Hcr \ min problem for all values A g] — oo, oo[. 

4.2. Separation of Variables for Partition of the State Space 

Another important property of separation of variables is valid for all divergences which have the form 
of a sum of convex functions f{pi,p*). Let the set of states be divided into two subsets, Ji and I2, 
and let the functionals n^, . . . be linear. We represent each probability distribution as a direct sum 
P = ® P^, where P^'^ are restrictions of P on Ji 2. 

Let us consider the problem 

H{P\\P*) min 

subject to conditions n*(P) = Ui for a set of linear functionals u^{P). 

The solution P™"' to this problem has a form P™'" = P^™'" ® P^^'"^, where P^'^ are solutions to 
the problems 

^(pi,2||p*i,2^) ^ j^in 

subject to conditions m'(Pi 2) = U^'"^ and Xlie/ia^'i^ ~ ^1.2 some redistribution of the linear 
functionals values, Ui = Ul + Uf, and of the total probability, 1 = tti + 7r2 {tii^2 > 0) . 

Again, the solution to the divergence minimization problem is composed from solutions of the partial 
maximization problems. Let us call this property the separation of variables for incompatible events 
(because Ji n /2 = 0). 

This property is trivially valid for the Tsallis family (for a > 0, and for a < with the change of 
minimization to maximization) and for the CR family. For the Renyi family it also holds (for a > 0, and 
for a < with the change from minimization to maximization), because the Renyi entropy is a function 
of those trace-form entropies, their level sets coincide. 
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Again, a simple check shows that this separation of variables property holds also for the convex 
combination of Shannon's and Burg's entropies, (5Dy^i^{P\\P*) + (1 - (5)Dv^i^{P*\\P). 

The question arises: is there any new divergence that has the following three properties: (i) the 
divergence H{P\\P*) should decrease in Markov processes with equilibrium P*, (ii) for minimization 
problems the separation of variables for independent subsystems holds and (iii) the separation of 
variables for incompatible events holds. A new divergence means here that it is not a function of a 
divergence from the CR family or from the convex combination of the Shannon and the Burg entropies. 

The answer is: no, any divergence which has these three properties and is defined and differentiable 
for positive distributions is a function of for h{x) = or h{x) = /3x In x — (1 — /3) In x. If we relax 
the differentiability property, then we have to add to the CR family the CR analogue of min-entropy 



The limiting case for the CR family for A — )• — oo is less known but is also a continuous and piecewise 
differential Lyapunov function for the Master equation: 



Both properties of separation of variables are based on the specific additivity properties: additivity 
with respect to the composition of independent systems and additivity with respect to the partitioning 
of the space of states. Separation of variables can be considered as a weakened form of additivity: not 
the minimized function should be additive but there exists such a monotonic transformation of scale 
after which the function becomes additive (and different transformations may be needed for different 
additivity properties). 

4.3. "No More Entropies" Theorems 

The classical Shannon work included the characterization of entropy by its properties. This meant 
that the classical notion of entropy is natural, and no more entropies are expected. In the seminal work of 
Renyi, again the characterization of entropy by its properties was proved, and for this, extended family 
the no more entropies theorem was proved too. In this section, we prove the next no more entropies 
theorem, where two one-parametric families are selected as sensible: the CR family and the convex 
combination of Shannon's and Burg's entropies. They are two branches of solutions of the correspondent 
functional equation and intersect at two points: Shannon's entropy (A = 1 in the CR family) and Burg's 
entropy (A = 0). We consider entropies as equivalent if their level sets coincide. In that sense, the Renyi 
entropy and the Tsallis entropy (with a > 0) are equivalent to the CR entropy with a — 1 = A, A>— 1. 

Following Renyi, we consider entropies of incomplete distributions: pi > 0, J2iPi — 1- The 
divergence H{P\\P*) is a smooth function of a pair of positive generalized probability distributions 



The following 3 properties are required for characterization of the "natural" entropies. 

1 . To provide the separation of variables for incompatible events together with the symmetry property 
we assume that the divergence is separable, possibly, after a scaling transformation: there exists 




£1 _1 




P = (Pt), Pi > and P* = (p*), p* > 0,1 = 1,. . .n. 
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such a function of two variables f{p,p*) and a monotonic function of one variable that 
H{P\\P*) = fiPuP*)). This formula allows us to define H{P\\P*) for all n. 

2. H{P\\P*) is a Lyapunov function for the Kolmogorov equation (22) for any Markov chain with 
equilibrium P*. (One can call these functions the universal Lyapunov functions because they do 
not depend on the kinetic coefficients directly, but only on the equilibrium distribution P*.) 

3. To provide separation of variables for independent subsystems we assume that H{P\\P*) is 
additive (possibly after a scaling transformation): there exists such a function of one variable ^/;(x) 
that the function ilj{H(P\\P*)) is additive for the union of independent subsystems: if P = (pij), 

p^J = wv p^j = then m{p\\p*)) = miQWQ*)) + m{R\m)- 

Theorem 1. If a -smooth divergence H{P\\P*) satisfies the conditions 1-3 then, up to 
monotonic transformation, it is either the CR divergence Hqr x or a convex combination of the 
Botlzmann-Gibbs-Shannon and the Burg entropies, Hh{P\\P*) = /3-Dkl(^'||-P*) + (1 - /3)^kl(^1I^)- 

In a paper [46] this family was identified as the Tsallis relative entropy with some abuse of language, 
because in the Tsallis entropy the case with a < is usually excluded. 

First of all, let us prove that any function which satisfies the conditions 1 and 2 is a monotone function 
of a Csiszar-Morimoto function (6) for some convex smooth function h{x). This was mentioned in 2003 
by P. Gorban [46]. Recently, a similar statement was published by S. Amari (Theorem 1 in [57]). 

Lemma 1. If a Lyapunov function H{p)for the Markov chain is of the trace-form (H {p) = J2i fiVii Pi) ) 
and is universal, then f{p,p*) = P*h{^) + const(p*), where h{x) is a convex function of one variable. 

Proof. Let us consider a Markov chain with two states. For such a chain 

dpi ,/p2 Pl\ ^ f Pi P2\ dp2 

If if is a Lyapunov function then H < and the following inequality holds: 

f df{p2,P*2) _ df{pi,pl) \ fPl_P2 \ ^ Q 
V 9p2 dpi J \pl pl) - 

We can consider pi,j92 as independent variables from an open triangle D = {{pi,P2) \ Pi,2 > 0, pi + 
P2 < 1}. For this purpose, we can include the Markov with two states into a chain with three states and 

gsi = qi3 = 0. 

If for a continuous function of two variables tlj{x,y) in an open domain C an inequality 
{'ip{xi,yi) — ip{x2,y2))iyi — y2) < holds then this function does not depend on x in D. Indeed, 
let there exist such values xi,2 and y that 'ip{xi,y) ^ ip{x2,y), 'ip{xi,y) — 'ip{x2,y) = e > 0. We can 
find such 5 > that {xi,y + Ay) e D and \ij{xi,y + Ay) - < e/2 if |Ay| < 5. Hence, 

i>{xi, y + Ay) — ^p{x2, y) > e/2 > if | Ay| < 5. At the same time {^p{xi, y + Ay) — ip{x2, y))Ay < 0, 
hence, for a positive < Ay < 5 we have a contradiction. Therefore, the function ^^^q^ •* is a monotonic 
function of ^, hence, f{p, p*) = p*h(^) + const(p*), where his a convex function of one variable. □ 
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This lemma has important corollaries about many popular divergences H{P{t)\\P*) which are not 
Lyapunov functions of Markov chains. This means that there exist such distributions Pq and P* and a 
Markov chain with equilibrium distribution P* that due to the Kolmogorov equations 

dH{P{t)\\P*) 

if -P(O) = Pq- This Markov process increases divergence between the distributions P, P* (in a vicinity 
of Po) instead of making them closer. For example, 

Corollary 1. The following Bregman divergences [58] are not universal Lyapunov functions for Markov 
chains: 

• Squared Euclidean distance B{P\\P*) = Xljfe ~PiY' 

• The Itakura-Saito divergence [59] B{P\\P*) = (f^ " log - l)- ^ 

These divergences violate the requirement: due to the Markov process distributions always 
monotonically approach equilibrium. (Nevertheless, among the Bregman divergences there exists a 
universal Lyapunov function for Markov chains, the Kulback-Leibler divergence.) 

We place the proof of Theorem 1 in Appendix. 

Remark. If we relax the requirement of smoothness and consider in conditions of Theorem 1 just 
continuous functions, then we have to add to the answer the limit divergences, 

^^CRoo(P||P*)=max|p|-l; 
HcK -oo(P|m= max 1^1-1 

5. Markov Order 

5.1. Entropy: a Function or an Order? 

Theorem 1 gives us all of the divergences for which (i) the Markov chains monotonically approach 
their equilibrium, (ii) the level sets are the same as for a separable (sum over states) divergence and (iii) 
the level sets are the same as for a divergence which is additive with respect to union of independent 
subsystems. 

We operate with the level sets and their orders, compare where the divergence is larger (for 
monotonicity of the Markov chains evolution), but the values of entropy are not important by themselves. 
We are interested in the following order: P precedes Q with respect to the divergence H,,,{P\\P*) if 
there exists such a continuous curve P{t) (t E [0, 1]) that P(0) = P, P(l) = Q and the function 
H(t) = iJ...(P(t)||P*) monotonically decreases on the interval t E [0, 1]. This property is invariant with 
respect to a monotonic (increasing) transformation of the divergence. Such a transformation does not 
change the conditional minimizers or maximizers of the divergence. 

There exists one important property that is not invariant with respect to monotonic transformations. 
The increasing function F(H) of a convex function H{P) is not obligatorily a convex function. 
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Nevertheless, the sublevel sets given by inequalities H{P) < a coincide with the sublevel sets 
F{H{P)) < F(a). Hence, sublevel sets for F(H(P)) remain convex. 
The Jensen inequality 

H{ep + (1 - e)Q) < eH{p) + (i - e)HiQ) 

(9 e [0,1]) is not invariant with respect to monotonic transformations. Instead of them, there appears the 
max form analogue of the Jensen inequality: 

H{eP+{l-e)Q) <max{H{P),H{Q)} , Be [0,1] (31) 

This inequality is invariant with respect to monotonically increasing transformations and it is equivalent 
to convexity of sublevel sets. 

Proposition 1. All sublevel sets of a function H on a convex set V are convex if and only if for any two 
points P,Q E V and every 9 G [0, 1] the inequality (31) holds. □ 

It seems very natural to consider divergences as orders on distribution spaces, and discuss only 
properties which are invariant with respect to monotonic transformations. From this point of view, the 
CR family appears absolutely naturally from the additivity (ii) and the "sum over states" (iii) axioms, 
as well as the convex combination f3DKL{P[\P*) + (1 - f3)DKL{P*[\P) (« G [0, 1]), and in the above 
property context there are no other smooth divergences. 

5.2. Description of Markov Order 

The CR family and the convex combinations of Shannon's and Burg relative entropies are 
distinguished families of divergences. Apart from them there are many various "divergences", and even 
the Csiszar-Morimoto functions (6) do not include all used possibilities. Of course, most users prefer 
to have an unambiguous choice of entropy: it would be nice to have "the best entropy" for any class of 
problems. But from some point of view, ambiguity of the entropy choice is unavoidable. In this section 
we will explain why the choice of entropy is necessarily non unique and demonstrate that for many 
MaxEnt problems the natural solution is not a fixed distribution, but a well defined set of distributions. 

The most standard use of divergence in many application is as follows: 

1. On a given space of states an "equilibrium distribution" P* is given. If we deal with the probability 
distribution in real kinetic processes then it means that without any additional restriction the current 
distribution will relax to P*. In that sense, P* is the most disordered distribution. On the other 
hand, P* may be considered as the "most disordered" distribution with respect to some a priori 
information. 

2. We do not know the current distribution P, but we do know some linear functional, the moments 

uiP). 

3. We do not want to introduce any subjective arbitrariness in the estimation of P and define it as the 
"most disordered" distribution for given value u{P) = U and equilibrium P*. That is, we define 
P as solution to the problem: 

H„XP\\P*) min subject to u{P) = U (32) 
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Without the condition u{P) = U the solution should be simply P*. 

Now we have too many entropies and do not know what is the optimal choice of H and what should 
be the optimal estimate of P. In this case the proper question may be: which P could not be such 
an optimal estimate! We can answer the exclusion question. Let for a given P° the condition hold, 
u{P^) = U. If there exists a Markov process with equilibrium P* such that at point P° due to the 
Kolmogorov equation (22) 

dt ^ dt 

then P° cannot be the optimal estimate of the distribution P under condition u{P) = U. 

The motivation of this approach is simple: any Markov process with equilibrium P* increases disorder 
and brings the system "nearer" to the equilibrium P*. If at P° it is possible to move along the condition 
plane towards the more disordered distribution then P° cannot be considered as an extremely disordered 
distribution on this plane. On the other hand, we can consider P° as a possible extremely disordered 
distribution on the condition plane, if for any Markov process with equilibrium P* the solution of the 
Kolmogorov equation (22) P(t) with initial condition P(0) = P° has no points on the plane u{P) = U 
for t > 0. 

Markov process here is considered as a "randomization". Any set C of distributions can be divided in 
two parts: the distributions which retain in C after some non-trivial randomization and the distributions 
which leave C after any non-trivial randomization. The last are the maximally random elements of 
C: they cannot become more random and retain in C. Conditional minimizers of relative entropies 
Hh{P\\P*) in C are maximally random in that sense. 

There are too many functions Hh{P\\P*) for effective description of all their conditional minimizers. 
Nevertheless, we can describe the maximally random distributions directly, by analysis of Markov 
processes. 

To analyze these properties more precisely, we need some formal definitions. 

Definition 1. (Markov preorder). If for distributions P° and P^ there exists such a Markov process with 
equilibrium P* that for the solution of the Kolmogorov equation with P(0) = P° we have P(l) = P^ 
then we say that P^ and P^ are connected by the Markov preorder with equilibrium P* and use notation 

po y% P\ 

Definition 2. Markov order is the closed transitive closure of the Markov preorder. For the Markov 
order with equilibrium P* we use notation P° >- p. P^. 

For a given P* = (p*) and a distribution P = (pi) the set of all vectors v with coordinates 

where p* and qij > are connected by identity (23) is a closed convex cone. This is a cone of all 
possible time derivatives of the probability distribution at point P for Markov processes with equilibrium 
P* = (p*). For this cone, we use notation Q(p,p*) 

Definition 3. For each distribution P and a n-dimensional vector A we say that A <(p,p*) ?/ A G 
Q(p p.). This is the local Markov order. 
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Proposition 2. Q(p,p*) is a proper cone, i.e., it does not include any straight line. 

Proof. To prove this proposition its is sufficient to analyze the formula for entropy production (for 
example, in form (24)) and mention that for strictly convex h (for example, for traditional xlnx or 
(x — 1)^/2) dHh/dt = if and only if dP/dt = 0. If the cone Q{p,p.) includes both vectors x and —a; 
(x 7^ it means that there exist Markov chains with equilibrium P* and with opposite time derivatives 
at point P. Due to the positivity of entropy production (24) this is impossible. □ 

The connection between the local Markov order and the Markov order gives the following proposition, 
which immediately follows from definitions. 

Proposition 3. P° >- p* P^ if and only if there exists such a continuous almost everywhere differentiable 
curve P{t) in the simplex of probability distribution that -P(O) = P^, -P(l) = P^ and for all t G [0, 1], 
where P{t) is differentiable, 

^ e Q(P(,),p.) □ (33) 

For our purposes, the following estimate of the Markov order through the local Markov order 
is important. 

Proposition 4. If P^ yp, P^ then P^ >(po,p.) P\ i.e., P^ - P^ e Q(p,p*). 

This proposition follows from the characterization of the local order and detailed description of the 
cone Q(p{t),p*) (Theorem 2 below). 

Let us recall that a convex pointed cone is a convex envelope of its extreme rays. A ray with directing 
vector X is a set of points Ax (A > 0). We say that / is an extreme ray of Q if for any u G / and any 
x,|/ G Q, whenever u = (x + ?/)/2, we must have x,y G /. To characterize the extreme rays of the 
cones of the local Markov order Q(p,p*) we need a graph representation of the Markov chains. We use 
the notation Ai for states (vertices), and designate transition from state Ai to state Aj by an arrow (edge) 
Ai — )■ Aj. This transition has its transition intensity qji (the coefficient in the Kolmogorov equation 
(21)). 

Lemma 2. Any extreme ray of the cone Q(p,p.) corresponds to a Markov process which transition graph 
is a simple cycle 

■Ail ~^ -A'h ~^ ■ ■ ■ -^ik ~^ ^ii 

where k < n, all the indices ii, . . .ik are different, and transition intensities for a directing vector of 
such an extreme ray qi.^^ i. may be selected as l/p*.: 

%+i ^j = ^ (34) 
(here we use the standard convention that for a cycle qi^^^ = qi^ ij. 

Proof. First of all, let us mention that if for three vectors x,y,u G Q(p,p*) we have u = [x + y)/2 
then the set of transitions with non-zero intensities for corresponding Markov processes for x and y are 
included in this set for u (because negative intensities are impossible). Secondly, just by calculation 
of the free variables in the equations (23) (with additional condition) we find that the the amount of 
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non-zero intensities for a transition scheme which represents an extreme ray should be equal to the 
amount of states included in the transition scheme. Finally, there is only one scheme with k vertices, k 
edges and a positive equilibrium, a simple oriented cycle. □ 

Theorem 2. Any extreme ray of the cone Q(p,p*) corresponds to a Markov process whose transition 
graph is a simple cycle of the length 2: Ai ^ Aj. A transition intensities qij, qjifor a directing vector 
of such an extreme ray may be selected as 

(lij = -T ' l3i = zi; (35) 
Pj Pi 

Proof. Due to Lemma 2, it is sufficient to prove that for any distribution P the right hand side of the 
Kolmogorov equation (22) for a simple cycle with transition intensities (34) is a conic combination (the 
combination with non-negative real coefficients) of the right hand sides of this equation for simple cycles 
of the length 2 at the same point P. Let us prove this by induction. For the cycle length 2 it is trivially 
true. Let this hold for the cycle lengths 2, ... n — 1. For a cycle of length n, Ai^ — Ai^ Ai^ — )• Ai^, 

with transition intensities given by (34) the right hand side of the Kolmogorov equation is the vector v 
with coordinates 

Pi 3-1 Pij 



Pi 



(under the standard convention regarding cyclic order). Other coordinates of v are zeros. Let us find the 
minimal value of Pijp*^ and rearrange the indices by a cyclic permutation to put this minimum in the 
first place: 

• \Pij\ Pii 
mm < — - > = — 

The vector f is a sum of two vectors: a directing vector for the cycle Ai^ — > . . . Ai^ — > Ai^ of the length 
n — \ with transition intensities given by formula (34) (under the standard convention about the cyclic 
order for this cycle) and a vector 



Pin 



P, 



where v"^ is the directing vector for a cycle of length 2, Ai^ which can have only two 

non-zero coordinates: 

2 Pi2 Ph 2 

Pi2 Pil 

The coefficient in front of v"^ is positive because Pi^/p*^ is the minimal value of PijP*^- A case when 
Pil/ Pil = Pi2/Pi2 ^^^^ need special attention because it is equivalent to the shorter cycle Ai^ — )■ 
Ai^ — )■ Ai^ (Ai^ could be omitted). A conic combination of conic combinations is a conic 
combination again. □ 

It is quite surprising that the local Markov order and, hence, the Markov order also are generated 
by the reversible Markov chains which satisfy the detailed balance principle. We did not include any 
reversibility assumptions, and studied the general Markov chains. Nevertheless, for the study of orders, 
the system of cycles of length 2 all of which have the same equilibrium is sufficient. 
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5.3. Combinatorics of Local Markov Order 

Let us describe the local Markov order in more detail. First of all, we represent kinetics of the 
reversible Markov chains. For each pair Aj, Aj {i ^ j) we select an arbitrary order in the pair and write 
the correspondent cycle of the length 2 in the form Ai ^ Aj. For this cycle we introduce the directing 
vector 7*^ with coordinates 

-Sik + Sjk (36) 



11 
Ik 



where 5ik is the Kronecker delta. This vector has the zth coordinate —1, the jth coordinate 1 and other 
coordinates are zero. Vectors 7*^ are parallel to the edges of the standard simplex in i?". They are 
antisymmetric in their indexes: 7*^ = —7-'*. 

We can rewrite the Kolmogorov equation in the form 



dP 



w 



(37) 



pairs ij 



where i ^ j, each pair is included in the sum only once (in the preselected order of i, j) and 



Pi 

P* 



The coefficient r^j > satisfies the detailed balance principle: 



QjiPi = QijPj 



We use the three-value sign function: 



signx 



0, 
1, 



if X < 0; 
if X = 0; 
if X > 



(38) 



With this function we can rewrite Equation (37) again as follows: 



dP_ 

lit ~ ^ 



pairs ij, fji^O 





p^_Pl 


\Pi P*J 


p*i Pj 



(39) 



The non-zero coefficients r^j may be arbitrary positive numbers. Therefore, using Theorem 2, we 
immediately find that the cone of the local Markov order at point P is 



Q{p,p*) = cone <^ 7'^sign — 



Pi 
P*i 



Pj 
Pj 



Tri > 



(40) 



Pi 



For a 



where cone{} stands for the conic hull. 

The number sign ( ^ — ^ ) is 1, when ^ > ^, —1, when ^ < ^ and 0, when ^ 

\Pi P] J Pi P] Pi Pj Pi 

given P*, the standard simplex of distributions P is divided by planes = ^ into convex polyhedra 

Pi Pj 

where functions sign — %] are constant. In these polyhedra the cone of the local Markov order (40) 



Q(p,p*) is also constant. Let us call these polyhedra compartments. 
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Figure 1. Compartments C^, corresponding cones (the angles) and all tableaus a for 
the Markov chain with three states (the choice of equilibrium {p* = 1/3), does not affect 
combinatorics and topology of tableaus, compartments and cones). 




In Figure 1 we represent compartments and cones of the local Markov order for the Markov chains 
with three states, ^1,2,3- The reversible Markov chain consists of three reversible transitions Ai ^ A2 ^ 
As ^ Ai with corresponding directing vectors 7^^ = (—1, 1, 0)"*^; 7^^ = (0, —1, l)"""; 7^^ = (1, 0, —1)^. 
The topology of the partitioning of the standard simplex into compartments and the possible values of 
the cone Q(p,p*) do not depend on the position of the equilibrium distribution P*. 

Let us describe all possible compartments and the correspondent local Markov order cones. For 
every natural number k < n — 1 the A;-dimensional compartments are numerated by surjective functions 
a : {1, 2, . . . , n} — )■ {1, 2, . . . , A; + 1}. Such a function defines the partial ordering of quantities ^ inside 
the compartment: 

a{j) (41) 



P* P* P* P* ^' 

Let us use for the correspondent compartment notation and for the Local Markov order cone 
Qcr- Let ki be a number of elements in preimage of i (i = 1, . . . , /c): ki = \{j \ a{j) = It 



Pi 



Pi 
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is convenient to represent surjection a as, a. tableau with k rows and ki cells in the ith row filled by 
numbers from {1,2,. First of all, let us draw diagram, that is a finite collection of cells arranged 
in left-justified rows. The ith row has ki cells. A tableau is obtained by filling cells with numbers 
{1,2,..., n}. Preimages of i are located in the 2th row. The entries in each row are increasing. (This is 
convenient to avoid ambiguity of the representation of the surjection a by the diagram.) Let us use for 
tableaus the same notation as for the corresponding surjections. 

Let a tableau A have k rows. We say that a tableau B follows A (and use notation A ^ B)ii B has 
k — 1 rows and B can be produced from A by joining of two neighboring rows in A (with ordering the 
numbers in the joined row). For the transitive closure of the relation — we use notation 

Proposition 5. rdQ^ = U^^? ^ 

Here rdU stands for the "relative boundary" of a set U in the minimal linear manifold which includes 

U. 

The following Proposition characterizes the local order cone through the surjection a. It is sufficient 
to use in definition of (40) vectors (36) with i and j from the neighbor rows of the diagram (see 
Figure 1). 

Proposition 6. For a given surjection a compartment and cone Q„ have the following description: 

C. = {P\^ = ^ ioT a(t) = aU) and ^ > ^ for a(j) = a(t) + ll (42) 
I Pj Pi Pj J 

Q„ = cone{7'^' | a{j) = a{i) + 1} □ (43) 
Compartment is defined by equalities ^ = ^ where i,j belong to one row of the tableau a and 

Pi Pj 

inequalities ^ > ^ where j is situated in a row one step down from i in the tableau (cr(j) = a{i) + 1). 

Pi Pj 

Cone Qcr is a conic hull of Yli=i ^i^i+i vectors 7*-'. For these vectors, j is situated in a row one step 
down from i in the tableau. Extreme rays of are products of the positive real half-line on vectors 7*-' 
(43). 

Each compartment has the lateral faces and the base. We call the face a lateral face, if its closure 
includes the equilibrium P*. The base of the compartment belongs to a border of the standard simplex 
of probability distributions. 

To enumerate all the lateral faces of a A;-dimensional compartment Ccr of codimension s (in C^,) we 
have to take all subsets with s elements in {1,2,..., A;}. For any such a subset J the correspondent 
k — s-dimensional lateral face is given by additional equalities = for a{j) = a{i) + 1, i G J. 

Pi Pj 

Proposition 7. All k — s-dimensional lateral faces of a k-dimensional compartment are in bijective 
correspondence with the s-element subsets J C {1,2, . . . , k}. For each J the correspondent lateral face 
is given in by equations 

— = ^ forall 2 G J and aij) = aU) + 1 □ (44) 

P*i Pj 
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The 1 -dimensional lateral faces (extreme rays) of compartment C^j are given by selection of one 
number from {1,2,...,/;;} (this number is the complement of J). For this number r, the correspondent 
1 -dimensional face is a set parameterized by a positive number a g]1, a,.], = 1/ Ylia{i)<r Pi'- 

— = a, for a(i) < r ; — = b, for a(i) > r : 
a > 1 > 6 > 0, a ^ p* + 6 5^ = 1 

i,(T(j)>r 

The compartment is the interior of the /c-dimensional simplex with vertices P* and Vr (r = 
1,2, .. .k). The vertex Vr is the intersection of the correspondent extreme ray (45) with the border of 
the standard simplex of probability distributions: P = Vr if 

Pi = Pi'^r, for a{i) < r; Pi = for a{i) > r (46) 

The base of the compartment Ca is ak — 1-dimensional simplex with vertices tv (r = 1, 2, . . . k). 
It is necessary to stress that we use the reversible Markov chains for construction of the general 
Markov order due to Theorem 2. 

6. The "Most Random" and Conditionally Extreme Distributions 

6.1. Conditionally Extreme Distributions in Markov Order 

The Markov order can be used to reduce the uncertainty in the standard settings. Let the plane L 
of the known values of some moments be given: m*(P) = f/j on L. Assume also that the "maximally 
disordered" distribution (equilibrium) P* is known and we assume that the probability distribution is P* 
if there is no restrictions. Then, the standard way to evaluate P for given moment conditions m*(P) = f/j 
is known: just to minimize H J^P\\P*) under these conditions. For the Markov order we also can define 
the conditionally extreme points on L. 

Definition 4. Let L be an affine subspace of R", E„ be a standard simplex in R". A probability 
distribution P E L\~\T,n is a conditionally extreme point of the Markov order on L if 

(P + Q(p,p.))nL = {P} 

It is useful to compare this definition to the condition of the extremum of a differentiable function H 
on L: gradif ±L. 

First of all, it is obvious that in the case when all the moments m*(P) are just some of the values Pi, 
then there exists only one extreme point of the Markov order on L, and this point is, at the same time, the 
conditional minimum on L of all Csiszar-Morimoto functions Hh{P) (6) (see, for example. Figure 2). 
This situation is unstable, and for a small perturbation of L the set of extreme points of the Markov order 
on L includes the intersection of L with one of compartments (Figure 3a). For the Markov chains with 
three states, each point of this intersection is a conditional minimizer of one of the CR divergences (see 
Fig. 3a). Such a situation persists for all L in general positions (Figure 3b). The extreme points of the 
family (3 D}ii^{P\\P*) + (1 — /3)-Dkl(-P*||-P) form an interval which is strictly inside the interval of the 
extreme points of the Markov order on L. For higher dimensions of L fl S„ the Markov order on L also 
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Figure 2. If the moments are just some of pi then all points of conditionally minimal 
divergence are the same for all the main divergences and coincide with the unique 
conditionally extreme point of the Markov order (example for the Markov chain with three 
states, symmetric equilibrium {p* = 1/3)) and the moment plane p2 =const. 




includes the intersection of L with some compartments, however the conditional minimizers of the CR 
divergences form a curve there, and extreme points of the family /3Dkl(-P||-P*) + (1 — /3)-Dkl(-P*||-P) 
on L form another curve. These two curves intersect at two points (A = 0, —1), which correspond to the 
BGS and Burg relative entropies. 
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Figure 3. The set of conditionally extreme points of the Markov order on the moment 
plane in two general positions. For the main divergences the points of conditionally minimal 
divergence are distributed in this set. For several of the most important divergences these 
minimizers are pointed out. In this simple example each extreme point of the Markov order 
is at the same time a minimizer of one of the iJcR a (A g] — oo, +oo[) (examples for the 
Markov chain with three states, symmetric equilibrium {p* = 1/3)). 




a) 

L 




b) 

6. 2. How to Find the Most Random Distributions ? 

Let the plane L of the known values of some moments be given: (P) = J^j '^)Pj = Ui(i = 1, . . .m) 
on L. For a given divergence H(P\\P*) we are looking for a conditional minimizer P: 

H{P\\P*) min subject to u\P) = Ui{i = 1, . . . m) (47) 

We can assume that H(P\\P*) is convex. Moreover, usually it is one of the Csiszar-Morimoto 
functions (6). This is very convenient for numerical minimization because the matrix of second 
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derivatives is diagonal. Let us introduce the Lagrange multipliers /ij = 1, . . . m) and write the system 
of equations (/iq is the Lagrange multiplier for the total probability identity Xli = ^ • 



dH 

''J ..=1 

n 

^ = f/i ; (48) 

i=i 

n 

Here we have n + m + 1 equations forn + m + 1 unknown variables {pj, /ij, /io). 

Usually is a convex function with a diagonal matrix of second variables and the method of choice 
for solution of this equation (48) is the Newton method. On the / + 1st iteration to find P'+^ = + AP 
we have to solve the following system of linear equations 



dH 



— — — Ap, = /io + > /iiM . - — 

,=1 dPjdps p=pi dpj 



^m;.Ap, = 0; (49) 

n 

For a diagonal matrix of the second derivatives the first n equations can be explicitly resolved. If for the 
solution of this system (49) the positivity condition + Apj > does not hold (for some of j) then we 
should decrease the step, for example by multiplication AP := 6'AP, where 

< 9 < min 



pJ+Ap,<0 \Api\ 

For initial approximation we can take any positive normalized distribution which satisfies the 
conditions n''(P) = Ui (i = 1, . . . m). 

For the Markov orders the set of conditionally extreme distributions consists of intersections of L 
with compartments. 

Here we find this set for one moment condition of the form u{P) = J2j '^jPj = U . First of all, assume 
that U ^ U*, where U* = u{P*) = J2j UjP* (if U = U* then equilibrium is the single conditionally 
extreme distribution). In this case, the set of conditionally extreme distributions is the intersection of the 
condition hyperplane with the closure of one compartment and can be described by the following system 
of equations and inequalities (under standard requirements Pi > 0, J2iPi = 1 )• 

J2 UjPj = U ; 

^ (50) 

if u.{U - W) > u,{U - W) 



Pi Pj 



(hence, ^ = ^ifui = Uj). 
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To find this solution it is sufficient to study dynamics of u{P) due to equations (37) and to compare it 



and, therefore, P* — P G Q(p,p*)- Equations and inequalities (50) mean that the set of conditionally 
extreme distributions is the intersection of the condition hyperplane with the closure of compartment C. 
In C, numbers ^ have the same order on the real line as numbers Ui(U — U*) have, these two tuples of 
numbers correspond to the same tableau a and C = C^. 

For several linearly independent conditions there exists a condition plane L: 



Let us introduce the m-dimensional space T with coordinates u\ Operator u{P) = (n*(P)) maps the 
distribution space into T and the affine manifold L (51) maps into a point with coordinates = Ui. 

If P* E L then the problem is trivial and the only extreme distribution of the Markov order on L is 
P*. Let us assume that P* ^ L. 

For each distribution P E Lwe can study the possible direction of motions of projection distributions 
onto T due to the Markov processes. 

First of all, let us mention that if u{Y-') = then the transitions Ai ^ Aj move the distribution 
along L. Hence, for any conditionally extreme distribution P E L this transition A^ ^ Aj should be in 
equilibrium and the partial equilibrium condition holds: ^ = ^. 



Let us consider processes with u{Y^) ^ 0. If there exists a convex combination (39) of vectors 
M(7*-')sign ( — ^ ) {u{j^^) 7^ 0) that is equal to zero then P cannot be an extreme distribution of the 

\ Pi Pj / 

Markov order on L. 

These two conditions for vectors 7*^ with u(Y'') = and for the set of vectors with non-zero 
projection on the condition space define the extreme distributions of the Markov order on the condition 
plane L for several conditions. 

7. Generalized Canonical Distribution 

7. 1 . Reference Distributions for Main Divergences 

A system with equilibrium P* is given and expected values of some variables Uj(P) = Uj are known. 
We need to find a distribution P with these values Uj{P) = Uj and is "the closest" to the equilibrium 
distribution under this condition. 

This distribution parameterized through expectation values is often called the reference distribution 
or generalized canonical distribution. After Gibbs and Jaynes, the standard statement of this problem is 
an optimization problem: 



for appropriate divergence H{P\\P*). If the number of conditions is m then this optimization problem 
can be often transformed into m + 1 equations with m + 1 unknown Lagrange multipliers. 

In this section, we study the problem of the generalized canonical distributions for single condition 



with dynamics of u{P) due to a model system P = P* — P. This model system is also a Markov chain 




(51) 



H{P\\P*) min, Uj{P) = Uj 



u{P) = E 



n 



UiPi = U,U^ U*. 



i=l 
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For the Csiszar-Morimoto functions Hh{P\\P* 



h'('-!l] (52) 



dpi \p* 

We assume that the function h'{x) has an inverse function g: g{h'{x)) = x for any x g]0, oo[. The 
method of Lagrange multipliers gives for the generalized canonical distribution: 



= '"o -i \- ^TT- , h \ — \ = IJ.O + iMi, 2_^pi = l, 2_^p^u^ = U (53) 

opi opi opi Wi J 

As a result, we get the final expression for the distribution 

Vi = Pi 9(1^0 + UifJ-) 
and equations for Lagrange multipliers /xq and 

n n 

^P^aifJ'O + Uin) = 1, ^p*g[no + Uin)ui = U (54) 

1=1 i=l 

If the image of h'{x) is the whole real line (h'(\0, oo[) = R) then for any real number y the value g{y) > 
is defined and there exist no problems about positivity of Pi due to (54). 

For the BGS relative entropy h'(x) = Inx (we use the normalized h{x) = xlnx — {x — 1) (19)). 
Therefore, g{x) = exp x and for the generalized canonical distribution we get 

n n n 

i=l 1=1 i=l 

As a result, we get one equation for ^ and an explicit expression for through p. 

These /iq and /i have the opposite sign comparing to (5) just because the formal difference between 
the entropy maximization and the relative entropy minimization. Equation (55) is essentially the same 
as (5). 

For the Burg entropy h'{x) = —^,g{x) = — ^ too and 

Pi = (56) 

/iO + Uifi 

For the Lagrange multipliers /io, /i we have a system of two algebraic equations 

n ^ n ^ 

y^^ = -i,y^^ = -u (57) 

^ /io + Uifl ~t 1^0 + Uill 

For the convex combination of the BGS and Burg entropies h\x) = /3 Inx — (0 < /3 < 1), and 
the function x = g{y) is a solution of a transcendent equation 

p\nx-^—^ = y (58) 

X 

Such a solution exists for all real y because this h'(x) is a (monotonic) bijection of ]0, 00 [ on the real 
line. 
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Solution to Equation (58) can be represented through a special function, the Lambert function [62]. 
This function is a solution to the transcendent equation 

we"" = z 

and is also known as W function, function or modified logarithm Im^ [36]. Below we use the main 
branch w = \mz for which Im^ > if 2 > and ImO = 0. Let us write (58) in the form 

lnx-- = -A (59) 

X 



where 5 = (1 - /3)//3, A = -y/i3. Then 

X = e "e 



A lm(<5eA) 



Another equivalent representation of the solution gives 

6 



X 



Indeed, let us take z = 5/x and calculate exponent of both sides of (59). After simple transformations, 
we obtain ze^ = 6e^. 

The identity Ima = In a — In Ima is convenient for algebraic operations with this function. Many 
other important properties are collected in [62] . 

The generalized canonical distribution for the convex combination of the BGS and Burg divergence 
is [36] 

_M_ (50) 

where Aj = — ^(/io + Uifi), 5 = (1 — /3)//3 and equations (54) hold for the Lagrange multipliers. 
For small 1 — /3 (small addition of the Burg entropy to the BGS entropy) we have 

P. = p: (e-^' + ^ - ^^e-^) + o{{l - P?) 

For the CR family ^x) = h\x) = ^-^±^, g{x) = (^^^g^)^ and 

, / A(A + l)(/io + n,^) + l V 
Pi = Pi /X , ,s (61) 



(A + 1) 

For A = 1 (a quadratic divergence) we easily get linear equations and explicit solutions for /iq and 
/i. If A = I then equations for the Lagrange multipliers (54) become quadratic and also allow explicit 
solution. The same is true for A = | and j but explicit solutions to the correspondent cubic or quartic 
equations are too cumbersome. 

We studied the generalized canonical distributions for one condition u{P) = U and main families 
of entropies. For the BGS entropy, the method of Lagrange multipliers gives one transcendent equation 
for the multiplier /ii and explicit expression for /iq as a function of /ii (55). In general, for functions 
Hh, the method gives a system of two equations (54). For the Burg entropy this is a system of 
algebraic equation (57). For a convex combination of the BGS and the Burg entropies the expression 
for generalized canonical distribution function includes the special Lambert function (60). For the CR 
family the generalized canonical distribution is presented by formula (61). for several values of A it can 
be represented in explicit form. The Tsallis entropy family is a subset of the CR family (up to constant 
multipliers). 
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7.2. Polyhedron of Generalized Canonical Distributions for the Markov Order 

The set of the most random distributions with respect to the Markov order under given condition 
consists of those distributions which may be achieved by randomization which has the given equilibrium 
distribution and does not violate the condition. 

In the previous section, this set was characterized for a single condition YliiPi'^i = U,U 7^ f/* by a 
system of inequalities and equations (50). It is a polyhedron that is an intersection of the closure of one 
compartment with the hyperplane of condition. Here we construct the dual description of this polyhedron 
as a convex envelope of the set of extreme points (vertices). 

The Krein-Milman theorem gives general backgrounds of such a representation of convex compact 
sets in locally convex topological vector spaces [63]: a compact convex set is the closed convex hull of 
its extreme points. (An extreme point of a convex set K is a point x E K which cannot be represented 
as an average x = \{y + z) iov y,z E K,y,z ^ x.) 

Let us assume that there are < n different numbers in the set of numbers Mj(f/—f/*) . Thereexists 
the unique surjection a : {1, 2, ... n} — {1, 2, ... A; + 1} with the following properties: a{i) < a{j) if 
andonlyifui(f/-f/*) > Uj{U -U*) (hence, a{i) = a(j) if and only if Ui(f/- f/*) = Uj{U -U*)). The 
polyhedron of generalized canonical distributions is the intersection of the condition plane J2i Pi^i = ^ 
with the closure of Cfj. 

This closure is a simplex with vertices P* and Vr (r = 1,2, .. .k) (46). The vertices of the intersection 
of this simplex with the condition hyperplane belong to edges of the simplex, hence we can easily find 
all of them: the edge [x, y] has nonempty intersection with the condition hyperplane if either u{x) > 

Uk.u{y) <U ox u{x) < Uk.u{y) > U. This intersection is a single point P if u{x) ^ u(y): 

p = Xx+{l-\)y, \= ^^y"^, (62) 

If u{x) = u{y) then the intersection is the whole edge, and the vertices are x and y. 

For example, if U is sufficiently close to U* then the intersection is a simplex with k vertices Wr 
(r = 1,2, .. .k). Each Wj. is the intersection of the edge \P*, Vr] with the condition hyperplane. 

Let us find these vertices explicitly. We have a system of two equations 



(63) 



*i i,cr{i)>r 

a Y Uip* + b ^ Uip* = U 

Position of the vertex Wr on the edge [P*, Vr] is given by the following expressions 

^ = a, for a{i) < r ; ^ = b, for cr(i) > r 

2^i,<j(i)>rPi l^i,a(i)<r^iPi ~^ l^i,a(i)<rPi l^i, a(i)>r ^iPi 



b= 1 



Yli,(7{i)>rPi ^i,a(i)<r'^iPi ~ Yli,a(i)<rPi Yli, a{i)>r '^iPi 
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If 6 > for all r then the polyhedron of generalized canonical distributions is a simplex with vertices 
Wr- If the solution becomes negative for some r then the set of vertices changes qualitatively and some 
of them belong to the base of Ca- For example, in Figure 3a the interval of the generalized canonical 
distribution (ID polyhedron) has vertices of two types: one belongs to the lateral face, another is situated 
on the basement of the compartment. In Figure 3b both vertices belong to the lateral faces. 

Vertices Wr on the edges [P*, Vr] have very special structure: the ratio Pi/p* can take for them only 
two values, it is either a or h. 

Another form for representation of vertices Wr (64) can be found as follows. Wr belongs to the edge 
[P*, Vr], hence, Wj. = XP* + (1 — \)vr for some A G [0, 1]. Equation for the value of A follows from the 
condition u{wr) = U : \U* + (1 — A)n(tv) = U. Hence, we can use (62) with x = P*, y = v^- 

For sufficiently large value oiU — U* for some of these vertices h loses positivity, and instead of them 
the vertices on edges [vr-, Wg] (46) appear. 

There exists a vertex on the edge [urji'g] if either ^(f^.) > Uk,u{vq) < U ov u{vr) < Uk.u{vq) > U. 
If u{vr) 7^ u{vg) then his vertex has the form P = Xvj. + (1 — \)vq and for A the condition u{P) = U 
gives (62) with x = Vr, y = Vq. If u{vr) = u{vq) then the edge [u{vr).,u{vq)] belongs to the condition 
plane and the extreme distributions are u{vr) u{vq). 

For each of Vr the ratio Pi/p* can take only two values: or 0. Without loss of generality we can 
assume that q > r. For a convex combination A^v + (1 — ^)vq (1 > A > 0) the ratio Pi/p* can take three 
values: Aa,. + (1 — \)aq (for a{i) < r), (1 — X)aq (for r < a{i) < q) and (for a{i) > q). 

The case when a vertex is one of the is also possible. In this case, there are two possible values of 
Pi/p*, it is either or 0. 

All the generalized canonical distributions from the polyhedron are convex combinations of its 
extreme points (vertices). If the set of vertices is {wr}, then for any generalized canonical distributions 
P = ^ XiWi (\i > 0, Aj = 1). The vertices can be found explicitly. Explicit formulas for the extreme 
generalized canonical distributions are given in this section: (64) and various applications of (62). These 
formulas are based on the description of compartment C^r given in Proposition 7 and Equation (46). 

8. History of the Markov Order 

8.1. Continuous Time Kinetics 

We have to discuss the history of the Markov order in the wider context of orders, with respect to 
which the solutions of kinetic equations change monotonically in time. The Markov order is a nice and 
constructive example of such an order and at the same time the prototype of all of them (similarly the 
Master Equation is a simple example of kinetic equations and, at the same time, the prototype of all 
kinetic equations). 

The idea of orders and attainable domains (the lower cones of these orders) in phase space was 
developed in many applications: from biological kinetics to chemical kinetics and engineering. A kinetic 
model includes information of various levels of detail and of variable reliability. Several types of building 
block are used to construct a kinetic model. The system of these building blocks can be described, for 
example, as follows: 
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1. The list of components (in chemical kinetics) or populations (in mathematical ecology) or states 
(for general Markov chains); 

2. The list of elementary processes (the reaction mechanism, the graph of trophic interactions or 
the transition graph), which is often supplemented by the lines or surfaces of partial equilibria of 
elementary processes; 

3. The reaction rates and kinetic constants. 

We believe that the lower level information is more accurate and reliable: we know the list of 
component better than the mechanism of transitions, and our knowledge of equilibrium surfaces is better 
than the information about exact values of kinetic constants. 

It is attractive to use the more reliable lower level information for qualitative and quantitative study 
of kinetics. Perhaps, the first example of such a analysis was performed in biological kinetics. 

In 1936, A.N. Kolmogorov [64] studied the dynamics of a pair of interacting populations of prey (x) 
and predator {y) in general form: 

x = xS{x,y), y = yW{x,y) 

under monotonicity conditions: dS{x,y) /dy < 0, dW{x,y) /dy < 0. The zero isoclines, the lines at 
which the rate of change for one population is zero (given by equations S{x, y) = or W{x, y) = 0), 
are graphs of two functions y{x). These isoclines divide the phase space into compartments (generically 
with curvilinear borders). In every compartment the angle of possible directions of motion is given 
(compare to Figure 1). 

Analysis of motion in these angles gives information about dynamics without an exact knowledge of 
the kinetic constants. The geometry of the zero isoclines intersection together with some monotonicity 
conditions give important information about the system dynamics [64] without exact knowledge of the 
right hand sides of the kinetic equations. 

This approach to population dynamics was further developed by many authors and applied to various 
problems [65,66]. The impact of this work on population dynamics was analyzed by K. Sigmund in 
review [67]. 

It seems very attractive to use an attainable region instead of the single trajectory in situations with 
incomplete information or with information with different levels of reliability. Such situations are typical 
in many areas of engineering. In 1964, F. Horn proposed to analyze the attainable regions for chemical 
reactors [68]. This approach was applied both to linear and nonlinear kinetic equations and became 
popular in chemical engineering. It was applied to the optimization of steady flow reactors [69], to 
batch reactor optimization by use of tendency models without knowledge of detailed kinetics [70] and 
for optimization of the reactor structure [71]. Analysis of attainable regions is recognized as a special 
geometric approach to reactor optimization [72] and as a crucially important part of the new paradigm of 
chemical engineering [73]. Plenty of particular applications was developed: from polymerization [74] 
to particle breakage in a ball mill [75]. Mathematical methods for study of attainable regions vary from 
the Pontryagin's maximum principle [76] to linear programming [77], the Shrink- Wrap algorithm [78] 
and convex analysis. 
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Figure 4. Attainable regions from an initial distribution for a linear system with three 
components Ai,A2, in coordinates Ci, C2 (concentrations of Ai, A2) (C3 = const — Ci — 
C2) [60]: for a full mechanism Ai ^ A2 ^ A3 ^ A\ (outlined region), for a two-step 
mechanism A\ ^ A2, Ai ^ A^ (horizontally shaded region) and for a two-step mechanism 
Ai ^ A2, A2 ^ A3 (vertically shaded region). Equilibrium is a*. The dashed lines are 
partial equilibria. 




The connection between attainable regions, thermodynamics and stoichiometric reaction mechanisms 
was studied by A.N. Gorban in the 1970s. In 1979, he demonstrated how to utilize the knowledge about 
partial equilibria of elementary processes to construct the attainable regions [60]. 

He noticed that the set (a cone) of possible direction for kinetics is defined by thermodynamics and 
the reaction mechanism (the system of the stoichiometric equation of elementary reactions). 

Thermodynamic data are more robust than the reaction mechanism and the reaction rates are 
known with lower accuracy than the stoichiometry of elementary reactions. Hence, there are two 
types of attainable regions. The first is the thermodynamic one, which use the linear restrictions and 
the thermodynamic functions [79]. The second is generated by thermodynamics and stoichiometric 
equations of elementary steps (but without reaction rates) [60,80]. 

It was demonstrated that the attainable regions significantly depend on the transition mechanism 
(Figure 4) and it is possible to use them for the mechanisms discrimination [81]. 

Already simple examples demonstrate that the sets of distributions which are accessible from a given 
initial distribution by Markov processes with equilibrium are, in general, non-convex polytopes [60,82] 
(see, for example, the outlined region in Figure 4, or, for particular graphs of transitions, any of the 
shaded regions there). This non-convexity makes the analysis of attainability for continuous time Markov 
processes more difficult (and also more intriguing). 
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This approach was developed for all thermodynamic potentials and for open systems as well [34]. 
Partially, the results are summarized in [14,61]. 

This approach was rediscovered by F.J. Krambeck [83] for linear systems, that is, for Markov chains, 
and by R. Shinnar and other authors [84] for more general nonlinear kinetics. There was even an open 
discussion about priority [86]. Now this geometric approach is applied to various chemical and industrial 
processes. 

8.2. Discrete Time Kinetics 

In our paper we deal mostly with continuous time Markov chains. For the discrete time Markov 
chains, the attainable regions have two important properties: they are convex and symmetric with respect 
to permutations of states. Because of this symmetry and convexity, the discrete time Markov order is 
characterized in detail. As far as we can go in history, this work was begun in early 1970s by A. Uhlmann 
and P.M. Alberti. The results of the first 10 years of this work were summarized in monograph [87]. A 
more recent bibliography (more than 100 references) is collected in review [88]. 

This series of work was concentrated mostly on processes with uniform equilibrium (doubly 
stochastic maps). The relative majorization, which we also use in Section 5, and the Markov order with 
respect to a non-uniform equilibrium was introduced by P. Harremoes in 2004 [89]. He used formalism 
based on the Lorenz diagrams. 

9. Conclusion 

Is playing with non-classical entropies and divergences just an extension to the fitting possibilities 
(no sense — just fitting)? We are sure now that this is not the case: two one-parametric families of 
non-classical divergences are distinguished by the very natural properties: 

1. They are Lyapunov functions for all Markov chains; 

2. They become additive with respect to the joining of independent systems after a monotone 
transformation of scale; 

3. They become additive with respect to a partitioning of the state space after a monotone 
transformation of scale. 

Two families of smooth divergences (for positive distributions) satisfy these requirements: the 
Cressie-Read family [29,30] 



If we relax the differentiability property, then we have to add to the the CR family two limiting cases: 




A g] — oo, oo[ 



and the convex combination of the Burg and Shannon relative entropies [34,35]: 





1 ; 
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HcK _oo(P|m =max|^l -1 

Beyond these two distinguished one-parametric families there is the whole world of the 
Csiszar-Morimoto Lyapunov functionals for the Master equation (6). These functions monotonically 
decrease along any solution of the Master equation. The set of all these functions can be used to reduce 
the uncertainty by conditional minimization: for each h we could find a conditional minimizer of Hh{p). 

Most users prefer to have an unambiguous choice of entropy: it would be nice to have "the best 
entropy" for any class of problems. But from a certain point of view, ambiguity of the entropy choice 
is unavoidable, and the choice of all conditional optimizers instead of a particular one is a possible 
way to avoid an arbitrary choice. The set of these minimizers evaluates the possible position of a 
"maximally random" probability distribution. For many MaxEnt problems the natural solution is not 
a fixed distribution, but a well defined set of distributions. 

The task to minimize functions Hh{p) which depend on a functional parameter h seems too 
complicated. The Markov order gives us another way for the evaluation of the set of possible "maximally 
random" probability distribution, and this evaluation is, in some sense, the best one. We defined the 
Markov order, studied its properties and demonstrated how it can be used to reduce uncertainty. 

It is quite surprising that the Markov order is generated by the reversible Markov chains which 
satisfy the detailed balance principle. We did not include any reversibility assumptions and studied 
the general Markov chains. There remain some questions about the structure and full description of the 
global Markov order. Nevertheless, to find the set of conditionally extreme ("most random") probability 
distributions, we need the local Markov order only. This local order is fully described in Section 5.2 
and has a very clear geometric structure. For a given equilibrium distribution P*, the simplex of 
probability distributions is divided by n{n — l)/2 hyperplanes of "partial equilibria" (this terminology 
comes from chemical kinetics [60,61]): = (there is one hyperplane for each pair of states 
In each compartment a cone of all possible time derivatives of the probability distribution is defined as a 
conic envelope of n(n — l)/2 vectors (40). The extreme rays of this cone are explicitly described in 
Proposition 6 (43). This cone defines the local Markov order. When we look for conditionally 
extreme distributions, this cone plays the same role as a hyperplane given by entropy growth condition 
( dS/ dt > 0) in the standard approach. 

For the problem of the generalized canonical (or reference) distribution the Markov order gives a 
polyhedron of the extremely disordered distributions. The vertices of that polyhedron can be computed 
explicitly. 

The construction of efficient algorithms for numerical calculation of conditionally extreme compacts 
in high dimensions is a challenging task for our future work as well as the application of this methodology 
to real life problems. 
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Appendix 

Proof of Theorem 1. The problem is to find all such universal and trace-form Lyapunov functions H for 
Markov chains, that there exists a monotonous function F, such that F(H(P)) = F{H{Q)) + F{H{R)) 

if P = Pij = QiTj. 

With Lemma 1 we get that 

Let F{x) and h(x) be differentiable as many times as needed. Differentiating the equality F{H{P)) = 
F{H{Q)) + F{H{R)) on ri and qi taking into account that g„ = 1 — Yl^=i 'ii f^m = ^ — '^T=^ 
we get that F'{H{P))H'^^^^{P) = -F"{H{P))H'JP)H^^{P), or, if = G{H{P)) then 

H' (P)H' (P) 

G{H{P)) = " ^ (65) 

qiri \ ) 
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It is possible if and only if every linear differential operator of the first order, which annulates H{P) and 
^Pi, annulates also 

f/;(p)f/;(p) 



H" (P) 

and it means that every differential operator which has the form 

/ dHjP) dH{P) \ d f dH{P) dH{P) \ d f dH{P) dH{P) \ _d 
V dq^ dqa J dq^ \ dq^ dq^ ) dq^ \ dq^ dqp ) dq 

= 3, 7 = 4 we get the following equation 



D 



i-f dqa 
annulates (66). For j3 = 2,a 



(66) 



-77- (67) 



F2{Q,R) [h' (fl 
F3{Q,R) 



h' r¥i 

94^1 



-h' 
-h' 
h' 



+ ^h" ( m 

_|_ q4rj_fi" I 94 ri 
qtr* \qlr* 



q2rm j^ll 
12^ m 



■h" 



qlr^ 



■h" 



92 '•m 
92 ^n 

qsr-m 

9|''*n 

qlrf, 



+ 



(68) 



where 



F,{Q,R) = Y,r, 

j 

F,{Q,R) = Y,r, 



qlr* 
qlr* 



h' 
h' 



F,{Q,R) = '£r,\h' 



h' 



qlr* 

q^ 

qlr* 

(Ml 
qlr* 



If we apply the differential operator ^ — ^> which annulates the conservation law = 1, to 

(68), and denote f{x) -- 
^, we get the equation 



the left part of (68), and denote fix) = xh"{x) + h'(x), Xi = ^, X2 = ^, x^ = ^,yi 

92 93 94 



'l ' m 



^2 



ifixsV-i) - f{x2yci) - fixsyi) + f{x2yi)){f{xiyi) - /(aJi^/a)) + 
ifixm) - fix^ys) - f\xiyi) + fixsy4)){f{x2yi) - f{x2y2)) + (69) 
{f{x2y-i) - /(xiyg) - f{x2yA) + f{xiy4)){f{x3yi) - /(XgJ/s)) = 
or, after differentiation on yi and y^ and denotation g(x) = f'{x) 

xig{xiyi){x3g{x-iy-i) - X2g{x2yz)) + X2g{x2yi){xig{xiy3) - (70) 
-X'igi^xzy'i)) + X35((x3?/i)(x25'(a;2?/3) - xigi^xiy^i)) = 
If 2/3 = 1, yi 7^ 0, (p{x) = xg{x), we get after multiplication (70) on yi 

(p{xiyi){(p{xs) - (p{x2)) + (p{x2yi)i(pixi) - (pix^)) + (p{xsyi){(p{x2) - ^(xi)) = (71) 

It implies that for every three positive numbers a, (3, 7 the functions Lp{ax), ip{(3x), fi'jx) are linearly 
dependent, and for Lp{x) the differential equation 

ax'^ip"{x) + bxip'{x) + cip{x) = (72) 

holds. This differential equation has solutions of two kinds: 
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1. ip{x) = Cix^'i + C2X^'^, ki 7^ k2, ki and k2 are real or complex-conjugate numbers. 

2. = Cix'^ + C2x'^ In X. 

Let us check, which of these solutions satisfy the functional equation (71). 

1. Lp{x) = Cix^^ + C2X^'2. After substitution of this into (71) and calculations we get 

CiC2{yi' - y^)(xt^x^^ - x^^x^^ + x^x^^ - x^x^^ + x^^x^^ - x^^x^^) = 

This means that Ci = 0, or C2 = 0, or ki = 0, or A;2 = and the solution of this kind can have 
only the form ip{x) = Cix'' + C2. 

2. Lp{x) = Cix^ + C2x'^ In X. After substitution of this into (71) and some calculations if yi 7^ we 
get 

Cl{{x\ - X2)x3lnX3 + (X3 - Xi)x2lnX2 + (X2 - X3)Xi Inxi) = 

This means that either C2 = and the solution is '^{x) = Cix^ or = and the solution is 

(y9(x) = Ci + C2 Inx. 

So, the equation (71) has two kinds of solutions: 

1. (^(x) = Cix'^ + C2, 

2. ip{x) = Ci + C2lnx 

Let us solve the equation /(x) = xh"{x) + h'{x) for each of these two cases. 

1. Lp{x) = Cix^ + C2, g{x) = Cix^~^ + there are two possibilities: 

1.1) A; = 0. Then ^(x) = f , /(x) = C\nx + Ci,h{x) = Cixlnx + C2lnx + C3X + C4; 

1.2) k 0. Then /(x) = Cx^ + Ci In x + C2, and here are also two possibilities: 

1.2.1) k = -1. Then/i(x) = Ci In^ x + Caxlnx + C3 Inx + C4X + C5; 

1.2.2) k ^ -1. Then/i(x) = Cix^+i + C2xlnx + Cglnx + C4X + C5; 

2. if{x) = Ci + C2lnx; g{x) = Ci^-^ + f ; /(x) = Ci In^ x + C2 Inx + C3; h{x) = Cixln^x + 
C2X In X + C3 In X + C4X + C5. 

(We have renamed constants during the calculations). 

For the next step let us check, which of these solutions remains a solution to equation (68). The result 
is that there are just two families of functions h{x) such, that equation (68) holds: 

1. h{x) = Cx^ + Cix + C2,k^0,k^l, 

2. h{x) = Cix In X + C2 In X + C3X + C4. 

The function h{x) should be convex. This condition determines the signs of coefficients Cj. 

The corresponding divergence H{P\\P*) is either one of the CR entropies or a convex combination 
of Shannon's and Burg's entropies up to a monotonic transformation. □ 

Characterization of Additive Trace-form Lyapunov Functions for Markov Chains. We will consider 
three important properties of Lyapunov functions H{P\\P*): 
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1. Universality: H is a Lyapunov function for Markov chains (22) with a given equilibrium P* for 
every possible values of kinetic coefficients > 0. 

2. if is a trace-form function. 

H{P\\P*) = Y,f{p-.P:) (73) 

i 

where / is a differentiable function of two variables. 

3. H is additive for composition of independent subsystems. It means that if P = pij = qiVj and 
P* = P*ij = q^r* then H{P\\P*) = H{Q\\Q*) + H{R\\R*). 

Here and further we suppose < Pi, p*, Qi, g*, rj, r* < 1. 

We consider the additivity condition as a functional equation and solve it. The following 
theorem describes all Lyapunov functions for Markov chains, which have all three properties 1) - 3) 
simultaneously. 

Let f(p, p*) be a twice differentiable function of two variables. 
Theorem 3. If a function H{P\\P*) has all the properties l)-3) simultaneously, then 

fiP^Pl = Plh (^) , H{P\\P*) = Y^plh {£j (74) 

where 

/i(x) = Cilnx + Csxlnx, Ci < 0, C2 > (75) 

Proof. We follow here the P. Gorban proof [46]. Another proof of this theorem was proposed in [90]. 
Due to Lemma 1 let us take H{P\\P*) in the form (74). Let h be twice differentiable in the interval 
]0, +00 [. The additivity equation 

H{P\\P*) - H{Q\\Q*) - H{R\\R*) = (76) 

holds. Here (in (76)) 

n— 1 m — 1 

Qn = -^Qi, rm = 1 - '^rj, P = Pij = QiTj 
i=l j=l 

Hipwn = J2^:r*h (^) , HiQim = Y^q^h , H{R\m = Y,r*h 

Let us take the derivatives of this equation first on qi and then on ri. Then we get the equation (g{x) = 
h'{x)) 

9i^,) - - 9{^) + aim + 

^1 1 1 ^1 m m 

Let US denote x = -2^^, y = -Sj^, z = -^i^, and ih(x) = q(x) + xq'(x). It is obvious that if n and m are 
more than 2, then x, y and z are independent and can take any positive values. So, we get the functional 
equation: 

^(?) =^(?/)+^W-^(^) (77) 
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Let's denote C2 = — V'(l) and ipi{a) = ip{a) — ipi^) and take x = 1. We get then 

the Cauchy functional equation [91]. The solution of this equation in the class of measurable functions is 
^i{a) = Ci Ina, where Ci is constant. So we get ^/^(x) = Cilnx + C2 and g{x) + xg'{x) = Cilnx + C2. 
The solution is g{x) = ^ + Cilnx + C2 - Cu h{x) = /(^ + Ci In x + C2 - Ci)dx = C3 In x + 
Cixlnx + {C2 — 2Ci)x + C4, or, renaming constants, h{x) = Ci Inx + C2x\nx + C^x + C4. In the 
expression for h{x) there are two parasite constants C3 and C4 which occurs because the initial equation 
was differentiated twice. So, C3 = 0, C4 = and h{x) = Cilnx + C2xlnx. Because h is convex, we 
have Ci<0 and C2 > 0. □ 

So, any universal additive trace-form Lyapunov function for Markov chains is a convex combination 
of the BGS entropy and the Burg entropy. 
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